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Abstract 

The statistics of energy levels for a disordered conductor are considered in 
the critical energy window near the mobility edge. It is shown that, if critical 
wave functions are multifractal, the one-dimensional gas of levels on the energy 
axis is compressible, in the sense that the variance of the level number in an 
interval is ((5N) 2 ) ~ x(N) for (N) 1. The compressibility, x = v/^d, is 
given exactly in terms of the multifractal exponent rj = d — Di at the mobility 
edge in a d-dimensional system. 
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An exact determination of energy levels in complex quantum systems is neither possible, 
nor desirable. Research in this area has focused instead on a statistical description of the 
spectra of such systems, which include: compound nuclei, quantum systems whose classi- 
cal behaviour is chaotic, and mesoscopic disordered conductors. In many cases, statistical 
characteristics turn out to be universal, i.e. independent of most microscopic details of the 
system and thus common to different systems which belong to the same universality class. 
In the following we shall be concerned with the specific example of a disordered conductor 
at an Anderson transition, for which universal spectral properties are an aspect of critical 
behaviour at the transition. For weakly disordered conductors, where wavefunctions corre- 
sponding to different levels have spatial overlap, there exists repulsion between the levels 
close in energy. Typically such repelling levels are governed [1-3] by Wigner-Dyson statistics 
(WDS) [4] which describe correlations between the eigenvalues in random-matrix ensembles 
of certain symmetries [5]. At the opposite extreme of strong disorder, spatial overlap be- 
tween wavefunctions corresponding to different levels is absent, the levels are uncorrelated 
and hence governed by Poisson statistics (PS). 

Repulsion between neighbouring levels is only one of the statistical characteristics which 
are totally different in WDS and PS. A second is spectral rigidity, which characterises fluc- 
tuations in the level density on scales large compared to the mean spacing. Specifically, 
consider the variance, S 2 ((iV)) = (iV 2 ) — (iV) 2 , for level number fluctuations in an energy 
interval which on average contains a large number, (N), of levels. Spectral rigidity is absent 
for PS: E 2 ((iV)) = (N), the maximal possible value. On the other hand, level fluctuations 
are strongly suppressed for WDS: £ 2 ((iV)) ~ In (AT) [5]; the levels are almost rigid. It is 
useful to think of the energy levels as the coordinates of a one-dimensional gas on the energy 
axis. If this gas has compressibility x, one expects X 2 ((iV)) ~ x{N) for (N) ^> 1. More 
precisely, we define (for system size L) 

x= > im .in,™). (1) 

(JV)^ooL-*oo d(N) V ' 

Clearly, x = 1 for a system of uncorrelated compressible levels (PS), and % = for a system 



of rigid, imcompressible levels (WDS). 

The subject we consider here is the existence of an intermediate universal limit for the 
compressibility (1). Spectral statistics in many systems experience a crossover from WDS to 
PS as some parameter of the system is changed. A smooth crossover of this type takes place 
for spectral statistics of electrons in 2d or quasi Id disordered conductors with increasing 
disorder. A new statistical regime [6,7] arises for d > 2 at the critical value of disorder 
that corresponds to the Anderson metal-insulator transition. The value of the spectral 
compressibility at the transition was first considered in Ref. [8] and has been a subject of 
controversy [7-10]. Note, however, that all the previous discussion has been limited to the 
framework of one-parameter scaling. 

We will show here that the one-parameter scaling is insufficient for the description of 
level statistics in the vicinity of the Anderson transition. The main result of the present 
work is an exact expression for the spectral compressibility at a mobility edge: 



Here r] = d — D 2 is one of the exponents characterising properties of the critical eigenstates: 
the set of multifractal dimensions, D p , govern the scaling behaviour of the moments of 
inverse participation ratio, 



for integer p. It is only the fractal dimensionality D 2 that enters Eq. (2), because the 
compressibility depends only on the two-level correlation function (TLCF), which (as we 
show below) can be expressed solely via second-order wavefunction correlations. 

Now we turn to the derivation of the result, Eq. (2). First, we express x Yld - the form 
factor of the TLCF. Then we outline a non-perturbative approach to calculating the form 
factor, developed previously in Refs. [11,12], which relates the spectral and wavefunction 
correlations. Finally, in the critical regime near the Anderson tranisition we relate the form 
factor to multifractal properties of the wavefunctions. 
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We define the TLCF as 



R(s) = (p}- 2 {(p(s + sA)p(s)}}-l, 



(4) 



where p(e) = L' d Y.n^{ e ~ E n ) is the density of states per unit volume, L is the sample 
size, and E n are eigenvalues of the Hamiltonian for a system of free electrons in a random 
potential. The mean density of states, (p), varies only on the scale of the Fermi energy, e F . 
Thus the mean level spacing A = l/((p) L d ) is considered constant. 

The number variance in an energy window of width E = (N)A <C e F is given in terms of 
the TLCF by E 2 ((iV)) = I%)((N) - \s\)R(s)ds. Differentiating with respect to (N), and 
taking the limit (N) — > oo after the thermodynamic limit, L — > oo, one obtains [7] 



and t H = h/A is the Heisenberg time. It is, of course, only in the thermodynamic limit 
that one can expect universal behaviour. Taking this limit at fixed (N) leaves A and t H , 
respectively, as the sole energy and time scales. Note that, for a system at a mobility edge, 
the number of critical levels (those for which the localisation length is bigger than the system 
size) diverges in the thermodynamic limit, as is required if there is to be a possibility that 
X takes a non-trivial critical value. We comment in passing that a well-know sum rule, 
f*™ R(s) ds = 0, can be derived if one takes limits in the opposite order: (N) — > oo before 
the thermodynamic limit, L — > oo. The fact that the order of limits does not necessarily 
commute has been discussed extensively elsewhere [9,10,13]. 

To obtain x, we use a recently developed [11,12] non-perturbative approach to calculating 
K(t) for times t <^t H . In this approach, K(t) is given by 




(5) 



where we have introduced the spectral form factor of the TLCF: 




(6) 



Here p(t) = (|^(0, t)| 2 ) is the ensemble-averaged quantum return probability for a 
wavepacket, ^(r,t) = v^oSn ^n{^)^n{f)^~ lErit ^ n , that initially occupies a small volume 
Vo near the origin [14]. The initial condition implies that the summation here is limited 
to the number of levels Mq ~ L d /V with energies lying within the energy band of width 
E ~ 1/pVb- Using spatial and spectral homogeneity, one can show that 



Thus p(t) and ultimately K(t), Eq. (7), are related to the local wavefunction correlations. 

Before applying Eq. (7) to find x, let us outline its derivation. A central idea is to 
consider the dependence of energy levels on some external parameter, A. Then, Eq. (7) can 
be obtained in two ways [11,12]. The first is phenomenological: we consider [11] a random 
walk in the ensemble of impurity configurations parametrised by A. Such a Brownian-motion 
level dynamics in the fictitious time oc A 2 is similar to that originally introduced by Dyson 
[15] in the context of random matrix theory. These dynamics generate the level correlations 
given by Eq. (7). Alternatively [12], using the homogeneity of level correlations at different 
points of the parametric space, one can prove for t ^ the identity 



which relates pit) to higher-order correlations between energy levels and wavefunctions (we 
stress that in constrast to Eq. (8), the summation here is extended over all eigenvalues of 
the Hamiltonian). 

In the limit t <C tu it is possible to factorise the average on the r.h.s. of Eq. (9). Consider 
first the sum over /. In both the metallic and the critical regimes, wavefunctions correspond- 
ing to different energies have spatial overlap, and so c„ in+ ; is a slowly decreasing function 
of \l\. Thus the oscillating factor, e - l ( E n-E n+ i)t'/n^ ac ^ g ag a cu t_ ff f or the sum over I: the 

number of terms contributing is of order tujt! > tn/t ^> 1. As a result, in the limit t — > 0, 
the entire contribution comes from terms for which |/| ^> 1. Now suppose that for one of 
these terms (i.e. fixed I, with |/| ^> 1) we carry out the sum on m. The result will depend on 
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the conditional distribution of E n+m , given that there are levels at the energies E n and E n+ \. 
Moreover, for t 7^ 0, non-zero contributions arise only from energy regions in which E n+m 
has a non-uniform distribution due to correlations with E n and E n+ i. The energy scale for 
such correlations is set only by A, and so the contribution to Eq.(10) comes from the two 
regions: \E n+m — E n \ < A and \E n+m — E n+t \ < A. Both sets of correlations make equal 
contributions to the r.h.s. of Eq. (9), as can be seen by changing variable in the integral 
from t' to t" = t — t'. Taking this into account, we factorise the average in Eq. (9) as 

tp (t) = 2 (^2 e-^-^™^ J* dt' (^2 c n ,n + ie l(En - En+l)t ' lf ^ , (10) 

This is exact for tjtu — > 0, provided that c ntn +i does not rapidly decrease with \l\ [16]. 

The first average in Eq. (10) reduces to K(t) on using Eqs. (4) and (6). The second 
average in Eq. (10) differs from the return probability (8) by a remainder arising from 
the terms with I > J\f . When \E n — E n+ i\ > E the two wavefunctions are essentially 
uncorrelated, and c ntn+ i = L~ d . The remainder is therefore 

/ E c n , n+l e-^~ E ^A = ( I p(E)e-^) = 2nhp5 (t) , 

Wo / XJE>Eo 1 

where 5o(t) is a <5-like function which is zero for t > t = h/E . As to/tn goes to zero in 
the thermodynamic limit considered, the integral in Eq. (10) is equal to rrhp + J * + p(t') dt'. 
Thefore, the factorisation (10) of the expression (9) leads to the relation (7) between K(t) 
and p(t) which is exact in the limit t/t H — > 0. 

Now we apply Eq. (7) to find the compressibility at the mobility edge. This regime is 
characterised by anomalous diffusion. This means that the diffusion coefficient D becomes 
time (and scale) dependent. In the absence of multifractality, it has been shown [17] by 
extending one-parameter scaling arguments to treat time- or frequency-dependent transport 
that D[t) oc t~ 1+2 l d . Then (restoring dimensional factors) p{t) ~ (Dt)~ d / 2 ~ hp/t. The 
numerator in Eq. (7) is therefore constant while the denominator is divergent in the universal 
limit, L — > 00 with t/t H fixed. Thus, in this case, K{t) = for t <^ t H and the compressibility 
(5) would be zero at the mobility edge, as in the metal. 
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However, there is no a priori reason to extend one-parameter scaling in this way. In 
fact, wavefunction correlations at the mobility edge are characterised by a set of multifractal 
exponents [18,19]. The relevant fractal dimension forp(t) is D 2 [20], Eq. (3), as it is expressed 
via \^ n \ 2 |^m| 2 , according to Eq. (8). Multifractality reveals itself in the time-dependence of 
pit). Since for I <ti E /A the factor c ntn+ i ~ \E n — E n+ i\~ v ^ d , we have [20]: 

p(f)~ VT' d (hp/t) 1 -*'* , V = d-D 2 . 

On subsituting this into Eq. (7), we note that the first, constant term in the denominator 
is negligible in the limit L — > oo, with t/tn fixed. Thus we obtain K(t) = 2i]/d for t <C t H . 

Note that in this limit R(s) = for s ^> 1, in agreement with the diagrammatic results 
of Ref. [7]. There are small t-dependent corrections to the leading behaviour of K(t) that 
govern the large- s asymptotics of R(s). These corrections, which will be discussed elsewhere, 
do not affect the limiting value K(0): thus our main result, Eq. (2), follows immediately. 

This result can be checked directly for d = 2. In this case, there is no Anderson tran- 
sition but a crossover from the weak to strong localisation regime. It is known that, in 
the weak localisation regime, where the dimensionless conductance, g = 2n 2 pD ^> 1, the 
high moments of the inverse participation ratio exhibit behaviour characterstic of that for 
critical states in d > 2. For the orthogonal, unitary and symplectic ensembles labelled by 
P = 1,2 or 4, the corresponding set of multifractal dimensions has been recently found 
non-perturbatively [21] to be D p = 2 —p/j3go, in agreement with the earlier renormalization 
group results [18,22]. This leads to rj = 2/(3g , and thus the limiting value of compressibility 
(2) is x — l/2/3<7o- R the other hand, the same value of compressibility has been obtained 
in this regime with the help of direct diagrammatic calculations [9] (the corresponding Eq. 
(15) in Ref. [9] contains an extra factor of 2 due to explicit inclusion of spin degeneracy). 

Another possibility [23] to check Eq.(2) is provided by a Id system with random hopping 
integrals Uj whose variance decreases as a power of the distance \i — j\: ((Uj) 2 ) = A\i — j\ 2a . 
This system is critical at a = 1, and eigenstates are multifractal. Calculations similar to 
those described above for 2D systems show [23] that Eq. (2) is again exactly satisfied. 
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Finally, the relation (2) between the compressibility x an d the multifractality exponent 
can be checked with the help of independent numerical calculations of both these quantities. 
Let us discuss available numerical results. For the 3D Anderson model (which belongs to 
the orthogonal symmetry class), the multifractality exponent D 2 was found [24,25] to be 
D 2 = 1.7 ± 0.2. Using Eq. (2) we have x — 0-22 ± 0.03 which is in a good agreement with 
the numerical result x = 0.2 0.3 [8,26]. The exponent rj has been studied numerically also 
for the Anderson transitions in 2D systems belonging to the unitary (quantum-Hall regime) 
and the symplectic symmetry classes [25,27]. In both cases rj 0.4 turns out to be small. 
Thus we can predict that the level compessibility should be very small, x ~ 0-1) i n the 
vicinity of these transitions. This could explain that in direct numerical simulations in the 
quantum-Hall regime [28] no linear contribution to S 2 has been found. 

In summary, we have derived an explicit relationship between the compressibility of en- 
ergy levels and the multifractal statistics of eigenfunctions in a disordered system near the 
Anderson transition. A non-zero compressibility implies that the one-parameter scaling is 
inapplicable to level statistics; for a hypothetical system where the one-parameter scaling 
description of the conductance could be extended to level correlations, no spectral compress- 
ibility would arise at the transition. The relationship (2) holds exactly for 2D critical states 
and is in a reasonable agreement with existing numerical results. 
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